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Competitive limits on the weakly interacting massive particle (WIMP) spin-independent scattering 
cross section are currently being produced by 76 Ge detectors originally designed to search for neu- 
trinoless double beta decay, such as the Heidelberg-Moscow and IGEX experiments. In the absence 
of background subtraction, limits on the WIMP interaction cross section are set by calculating the 
upper confidence limit on the theoretical event rate, given the observed event rate. The standard 
analysis technique involves calculating the 90% upper confidence limit on the number of events in 
each bin, and excluding any set of parameters (WIMP mass and cross-section) which produces a 
theoretical event rate for any bin which exceeds the 90% upper confidence limit on the event rate for 
that bin. We show that, if there is more than one energy bin, this produces exclusion limits that are 
actually at a lower degree of confidence than 90%, and are hence erroneously tight. We formulate 
criteria which produce true 90% confidence exclusion limits in these circumstances, including calcu- 
lating the individual bin confidence limit for which the overall probability that no bins exceeds this 
confidence limit is 90% and calculating the 90% minimum confidence limit on the number of bins 
which exceed their individual bin 90% confidence limits. We then compare the limits on the WIMP 
cross-section produced by these criteria with those found using the standard technique, using data 
from the Heidelberg-Moscow and IGEX experiments. 

PACS numbers: 98.80Cq, 14.80Ly astro-ph/yymmddd 



I. INTRODUCTION 

The non-Keplerian rotation curves of spiral galaxies 
imply that, provided that Newtonian gravity is correct, 
they are surrounded by halos of dark matter [0. The 
nature of the dark matter is unknown with possi- 
ble candidates including MAssive Compact Halo Objects 
(MACHOs), such as brown dwarves, Jupiters or black 
holes, and elementary particles, known as Weakly Inter- 
acting Massive Particles (WIMPs). Currently the best 
motivated WIMP candidate is the neutralino, the light- 
est supersymmetric particle |3j. 

Neutralinos, and WIMPs generically, can potentially 
be directly detected via their elastic scattering off target 
nuclei. Over the past five years the DAMA collaboration 
has found an annual modulation in the signal in their 
Nal detector, consistent with nuclear recoils caused by 
WIMPs with mass m x ~ 54 GeV. Their claim has been 
met with some scepticism |5|, due to the difficulties in- 
volved in extracting the small (~ a few per-cent) annual 
modulation signal from the backgrounds present, along 
with various technical concerns. It is therefore crucial to 
verify this result using other detectors and techniques. 

Currently the most stringent spin-independent WIMP 
exclusion limits come from the CDMS experiment's 10.6 
kgd exposure at Stanford (|]. They use the measured 
nuclear recoils in Si detectors (which are more sensitive 
to neutrons than to WIMPs), along with the multiple 
scattering events in their Ge detectors, to subtract the 
neutron background from their Ge single scatter events. 



They find, using a sophisticated likelihood analysis along 
with Monte-Carlo simulations (see Ref. j?J for details), 
that the Ge nuclear recoil events observed are most likely 
due to neutrons rather than WIMPs. The CDMS experi- 
ment will be moving to the Soudan mine, where the lower 
backgrounds should, in the long term, allow the improve- 
ment of the constraints on the WIMP cross-section by 
several orders of magnitude || . 

Competitive constraints are also coming from 76 Ge de- 
tectors originally designed to search for neutrinoless dou- 
ble beta decay, such as the Heidelberg-Moscow Q and 
IGEX jl0| experiments. The Heidelberg-Moscow experi- 
ment operates five 76 Ge detectors with total active mass 
of 10.96 kg, and to date have analysed data taken over 
0.249 years, from one detector with an active mass of 
2.758kg and 86% enrichment in 76 Ge, giving a total ex- 
posure of 250.836 kgd. The IGEX collaboration have to 
date analysed data taken over 30 days by one of their 
enriched detectors with an active mass of 2.0 kg, giving 
a total exposure of 60.0 kgd. In this paper we exam- 
ine the method used to calculate exclusion limits from 
these experiments. We outline the calculation of upper 
confidence limits on the number of events per bin, and 
formulate criteria which produce true overall 90% con- 
fidence exclusion limits on the WIMP scattering cross- 
section. We then compare the exclusion limits on the 
cross-section calculated using these criteria with those 
calculated using the standard analysis employed by the 
Heidelberg-Moscow and IGEX collaborations. 
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II. CONFIDENCE LIMIT CALCULATION 

In experiments without background subtraction, such 
as Heidelberg-Moscow (HM) § and IGEX @, no con- 
straints on the WIMP parameters can be derived from 
the shape of the observed energy spectrum. However any 
set of WIMP parameters (mass and cross-section) which 
would produce more events than are observed, at some 
confidence level, can be excluded at that confidence level. 
The probability of finding an observed number of events 
iVobs given an expected number of events A is given by 
the Poisson distribution: 



P{N obs \X) 



N obB \ 



(1) 



where the '|' denotes 'given'. Note that whilst iVobs 1S > 
by definition, an integer, A can take any positive value. 
Using Bayes theorem 
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where P(A) is the prior probability distribution for A. 
Following the Particle Data Group we take P(A) to be 
uniform in A |^] and the above expression then simplifies 
to 



P(A|iV obs ) =P(7V obs |A). 



(3) 



It then follows that the probability that A is greater than 



some value A p is 



P(A> Ap|7V obs ) 



P(X\N ohs )dX 



■dA. 



Using the relation 12 
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we find that 



No 
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= P(JV< AUs|A p ). 



(6) 



The HM and IGEX collaboration bin their data in bins 
of width 1 keV (42 bins from 9 keV to 51 kcV for HM 



*See Ref. |nj and references therein for a contrary 
however. 



and 45 bins from 4 keV to 49 keV for IGEX). In their 
analysis |}|,[l0| , the 90% confidence exclusion limit on the 
cross-section, cr p , is calculated by finding for each m x 
the value of <7 P for which the theoretical event rate in 
one of the bins just exceeds the 90% bin confidence limit 
on the event rate in that bin. The Heidelberg-Moscow 
collaboration employ a slight variation of this technique, 
re-binning the data in overlapping 5keV bins (i.e. 9-14 
keV, 10-15 keV etc.), claiming that otherwise the limits 
obtained are too conservative ||. 

By definition a 90% upper confidence limit means that 
there is a 10% probability that a theoretical number of 
events greater than the 90% confidence limit produced 
the observed number of events, for any given bin. We 
will now consider an ensemble of Nt bins. The probability 
distibution of the number of bins for which the theoretical 
number of events exceeds the 90% bin confidence limit, 
N e , is given by 



P{N C > 0) = 1 - P{N C = 0) = 1 - (0.9) Ni 
>0.1 if JV t > 1 . 



(7) 



Therefore the exclusion limits found using the standard 
analysis actually correspond to a lower degree of confi- 
dence than 90% and are hence erroneously tight. The 
larger the total number of bins the larger the discrep- 
ancy will be. For an ensemble of bins the probability 
that the observed number of events in iV c of the bins 
was produced by a theoretical number of events greater 
than the 100p% bin upper confidence limit is given by 
the binomial distribution: 



P(N C ) = C$j> N ' (1 - P ) 



Nt-Ne 



where 



C 



N t 

N B 



N c \(N t -N e )\' 



(8) 



(9) 



provided that each bin is independent [j] (this is not the 
case for the HM experiment's overlapping 5keV bins). 

We will now use the probability distribution of the 
number of bins N c exceeding their 100p% bin confidence 
limit (eq. (|J)) to formulate criteria which produce true 
90% minimum confidence exclusion limits. Firstly, for 
any given total number of bins, we can find the confi- 
dence level cl for which the probability that none of the 
bins exceed their 100cl% bin confidence is 90%: 



C^(l-cl) 



N t 



0.1. 



(10) 



' Whilst the values of A predicted for each bin by an underly- 
ing theoretical model are not independent of each other, the 
number of events observed in each bin is an independent Pois- 
son process, governed by the underlying value of A for that 
bin. 
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We can also calculate the 90% minimum upper con- 
fidence limit on the number of bins which exceed their 
90% bin confidence limit. We define A/90 as the smallest 
integer which satisfies: 



ATtot 

E 

Ne = Ng: 



C$0.1 



(0.9) 



(Nt-N.) 



< 0.1 . 



(11) 



Any set of WIMP parameters which produces N c > Ago 
is excluded at at least 90% confidence. For N t = 42 or 
45, V 90 = 8 and P(N C < V 90 ) = 0.946 (0.924) for N t = 
42(45). Since N c can only take on integer values this 
criteria does not produce exact 90% exclusion limits, but 
the amount by which the confidence limits are stronger 
than 90% is known, and fixed for fixed N t . 

A simple minded way to avoid the problem of calcu- 
lating 90% overall confidence limits would be to discard 
all but one of the energy bins; N t = 1 then and the 90% 
bin confidence limit gives an overall 90% confidence limit. 
The obvious choice for which energy bin to use is the low- 
est, threshold, energy bin, since di?/di?R decreases expo- 
nentially with increasing for the standard Maxwellian 
halo model. It is likely that this wasteful method will 
produce weaker exclusion limits than can be found using 
the entire data set though. 



III. EXCLUSION LIMITS 

We will now calculate the exclusion limits resulting 
from the HM and IGEX data, assuming a standard 
Maxwellian halo, using a variety of criteria: the HM 
collaboration's overlapping 5keV bins technique (crite- 
ria A-HM only), no bins exceeding their 90% bin confi- 
dence limit (criteria B), no bins exceeding their 100cl% 
bin confidence limit (criteria C), less than Ngo = 8 bins 
exceeding their 90% bin confidence limit (criteria D) and 
from the lowest energy bin not exceeding its 90% bin con- 
fidence limit (criteria E). The bin confidence limit cl for 
criteria C has been calculated from the binomial distri- 
bution to give 90% overall (as opposed to individual bin) 
exclusion limits; for the Nt = 42 bins of the HM experi- 
ment cl = 0.9975 whilst for the N t = 45 bins of the IGEX 
experiment cl = 0.9977. In Section || we found that Cri- 
teria D produces 94.6% confidence exclusion limits for 
the HM experiment and 92.4% confidence exclusion lim- 
its for the IGEX experiment. Note that A p /Ar o bs, for any 
fixed p, decreases with increasing N ^ s ; the smaller the 
number of counts observed the greater the uncertainty 
in the underlying see e.g. Refs. |Jl!|||3. theoretical 
number of counts. For details of the calculation of the 
theoretical event rate The resulting exclusion limits are 
plotted in Fig. 1, and in Tables I and II we list the en- 
ergy bins from which the exclusion limits on a p arise, as 
a function of m x , for each criteria. 

The constraints from criteria A and B are similar and 
are tighter than those found using either of the crite- 
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FIG. 1. The spin-independent WIMP cross-section exclu- 
sion limits found from the HM (top panel) and IGEX (lower 
panel) data, assuming the standard Maxwellian halo model, 
from: using the sliding 5keV bins method (long dashed line), 
no bin exceeding its 90% bin confidence limit (solid), no bin 
exceeding its 99.75% bin confidence limit (dotted), no more 
than 7 bins exceeding their 90% bin confidence limit (short 
dashed) and the lowest energy bin not exceeding its 90% bin 
confidence limit (dot-dashed) 



ria (C and D) calculated from the binomial distribu- 
tion to give 90% overall confidence limits, by roughly 
A [log 10 (Ccr)] ~ 0.2. As shown above (Eq. (g)) the re- 
quirement that no bin exceeds its 90% bin confidence 
limit (i.e. criteria B) actually produces exclusion limits 
at 1.2% (0.9%) confidence for N t = 42 (45) bins. Criteria 
D produces significantly weaker constraints than C for 
m x < 80 GeV, and we will discuss the reason for this 
below. For m x > 80 GeV the exclusion limit from crite- 
ria D is slightly weaker than that from criteria C since 
criteria D produces stronger than 90% confidence exclu- 
sion limits (94.6% for HM and 92.4% for IGEX). For the 
HM data criteria E produces exclusion limits which are 
identical to those from Criteria B for small m x , where 
the lowest energy bin is the most constraining, for large 
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m x (GeV) 


Energy bin(s) 


Criteria A 


< 54 


9-14 keV 


54-72 


11-16 keV 


> 72 


15-20 keV 


Criteria B 


< 44 


9-10 keV 


44-65 


11-12 keV 


65-81 


15-16 keV 


> 81 


18-19 keV 


Criteria C 


< 51 


9-10 keV 


51-65 


11-12 keV 


65-83 


12-13 keV 


83-142 


15-16 keV 


> 142 


18-19 keV 


Criteria D 


< 47 


12 3 4 5 6 7 


47-72 


1 2 3 4 6 7 10 


72-83 


1 3 4 6 7 8 10 


83-93 


3 4 6 7 8 9 10 


> 93 


3 4 7 8 9 10 11 



TABLE I. The most constraining energy bins for the Hei- 
delberg-Moscow data, for the criteria described in the text, 
for a standard Maxwellian halo. 

m x it produces the weakest exclusion limit since the more 
constraining higher energy bins have been discarded. For 
the IGEX data it turns out that the lowest energy bin 
is never the most constraining, and even for small m x 
criteria E produces weaker limits than Criteria B. 

For both data sets the number of events per bin 
initially decreases fairly sharply with increasing energy 
(Ey < 18(11) keV for HM (IGEX)). For larger ener- 
gies the number of events per bin varies between 2-17 
(0-10) for HM (IGEX). For the standard Maxwellian 
halo dR/dEft decreases exponentially with increasing 
E-r. The decrease is steepest for small m x , with 
dR/dER(Eix — > 0) being largest for small m x . For small 
m x the theoretical energy spectrum falls off more rapidly 
than the observed energy spectrum and the exclusion 
limit comes from the lowest energy bins, where the ob- 
served number of events is largest. The larger the ob- 
served number of events the smaller the difference be- 
tween Ago and A c i, and hence the difference in the ex- 
clusion limits produced by criteria B and C is smallest 
for small m x . For a rapidly decreasing theoretical energy 
spectrum a p has to be so large for eight bins to exceed 
their bin confidence limit that the theoretical event rate 
in the lowest bin is far greater than A c i, and criteria D 
produces misleadingly weak exclusion limits. 

If the WIMP velocity distribution is not close to 
Maxwellian, then the form of di?/d£ , R will deviate from 
exponential and the error in the exclusion limits pro- 
duced by the standard analysis technique could be more 



m x (GeV) 


Energy bin(s) 


Criteria B 


< 41 


5-6 keV 


41-131 


8-9 keV 


131-186 


14-15 keV 


> 186 


24-25 keV 


Criteria C 


< 47 


5-6 keV 


> 47 


8-9 keV 


Criteria D 


< 42 


12 3 4 5 6 7 


42-56 


1 2 3 4 5 6 11 


56-57 


1 2 3 5 10 11 12 


57-63 


1 2 4 5 10 11 12 


63-79 


2 4 5 6 10 11 12 


79-114 


2 4 5 10 11 12 21 


114-187 


2 5 10 11 12 14 21 


> 187 


4 5 10 11 12 14 21 



TABLE II. The most constraining energy bins for the 
IGEX data, for the criteria described in the text, for a stan- 
dard Maxwellian halo. 

significant. As an example we will study the exclusion 
limits for Sikivie's late infall model |l4|], where the CDM 
distribution at the Earth's location consists of a number 
of velocity flows plus a smooth isothermal background 
distribution. The resulting WIMP differential event rate 
dR/dEft has a series of steps and falls off more rapidly 
than that produced by a pure Maxwellian halo Jl5[ . 
Whilst it is unlikely that such velocity flows are present 
in the Milky Way , the real velocity distribution pro- 
duced by the hierarchical formation of the galactic halo is 
likely to deviate significantly from a smooth Maxwellian 
distribution, resulting in a dR/dE^ which varies signifi- 
cantly from the exponential decline produced by the stan- 
dard Maxwellian halo. The parameters of the late infall 
model can be found in Ref . jl5| , the total local halo den- 
sity is p x = 0.57 GeV cm -3 which corresponds to £ = 1.9. 
In Fig. 2 we plot the exclusion limits found from the HM 
data for this model, using the criteria defined above. In 
table III we list the bins from which the limit on a p arises 
for each criteria, as a function of m x . 

For m x < 50 GeV the difference in the exclusion lim- 
its produced by the different criteria is relatively small. 
As m x is increased the step-like drops in dR/dEn caused 
by the velocity flows move to higher Er [ p"5| , producing a 
sudden change in the predicted number of events for some 
of the bins, changing the bin from which the constraint 
on (Tp originates and, for criteria B and C, producing a 
kink in the exclusion limit curve. Criteria B produces 
limits which are tighter than those from Criteria C by 
roughly A [log 10 (Cc)] ~ 0.2, with the difference decreas- 
ing with decreasing m x , once more. For criteria C the 
cusps in the exclusion limit are shifted to slightly higher 
m x , Am x ~ 5 GeV for m x > 80 GeV, since A c i/Ago in- 
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FIG. 2. The spin-independent WIMP cross-section exclu- 
sion limits found for Sikivie's late infall model from the HM 
data found from: the sliding 5keV bins method (long dashed 
line), no bin exceeding its 90% bin confidence limit (solid), 
no bin exceeding its 99.75% bin confidence limit (dotted), 
no more than 7 bins exceeding the 90% bin confidence limit 
(short dashed) and the lowest energy bin not exceeding its 
90% bin confidence limit (dot-dashed). 

creases with decreasing iV b s - Once more criteria D pro- 
duces weaker exclusion limits then criteria C for m x < 80 
GeV, although for this halo model the difference between 
criteria C and D at large m x , is larger than for the stan- 
dard halo model. The cusps in the exclusion limits are 
smoothed out to some extent by criteria A and D since 
for these criteria the limits arise from a range of energies 
(a 5keV bins for A and 7 IkeV bins for D), rather than 
a single IkeV bin. 



IV. CONCLUSIONS 

We have shown that the standard analysis technique 
for calculating exclusion limits from WIMP direct detec- 
tion experiments without background subtraction pro- 
duces exclusion limits that are actually at a lower degree 
of confidence than the stated 90%, and are hence erro- 
neously tight, if there is more than one energy bin. 

We have examined criteria which produce true 90% 
minimum confidence exclusion limits for an ensemble of 
bins. The first (criteria C above) involves calculating 
the individual bin confidence limit for which the overall 
probability that no bin exceeds this confidence limit is 
90%. For the 42 bins of the HM experiment this indi- 
vidual bin limit is 99.75% whilst for the 45 bins of the 
IGEX experiment it is 99.77%. The second (D above) 
involves calculating the minimum 90% upper confidence 
limit on the number of bins which exceed their individual 
bin 90% confidence limit, for 42 or 45 bins this limit is 
8. Finally the third (E above) involves using only the 
lowest, threshold, energy bin. 



m x (GeV) 


Energy bins 


Criteria A 


< 93 
93-130 
130-141 
> 141 


9-14 keV 
11-16 keV 

14- 19 keV 

15- 20 keV 


Criteria B 


< 65, 80-85 
65-78 85-106 
78-80 106-124 
124-160 
> 160 


9-10 keV 

11- 12 keV 

12- 13 keV 
15-16 keV 
18-19 keV 


Criteria C 


< 68, 73-89 
68-73, 89-106 
106-128 
128-167 
> 167 


9-10 keV 

11- 12 keV 

12- 13 keV 
15-16 keV 
18-19 keV 


Criteria D 


< 43, 46-47, 66-87 
43-46, 48-59, 87-96, 108-113 
47-48, 59-66 
96-108, 
113-116 
116-135 
135-163 
163-190 
> 190 


12 3 4 5 6 7 
1 2 3 4 6 7 10 
1 2 3 4 5 7 10 

1 2 3 4 7 8 10 

2 3 4 5 6 7 10 

3 4 5 6 7 8 10 

3 4 6 7 8 9 10 

4 6 7 8 9 10 11 
4 7 8 9 10 11 12 



TABLE III. The most constraining energy bins for the 
Heidelberg-Moscow data, for the criteria described in the text, 
for Sikivie's late infall halo model. 

We then compared the exclusion limits produced by 
these criteria with those produce by the standard analy- 
sis technique, for both the HM and IGEX experiments for 
a standard Maxwellian halo, and for the HM experiment 
for Sikivie's late infall model. We found that the standard 
technique produces limits which are erroneously tighter 
by roughly A [log 10 «cr)] ~ 0.2. for m x > 80 GeV. The 
difference would be larger for experiments with larger 
numbers of bins and/or less events in the most constrain- 
ing bins. Criteria D produces exclusion limits which are 
misleadingly weak for small m x and does not produce 
exactly 90% exclusion limits, since the number of bins is 
an integer. Criteria E is wasteful and produces substan- 
tially weaker exclusion limits than could be found using 
the entire data set. We therefore recommend that crite- 
ria C should be used to produce overall 90% confidence 
limits. 

For halo models with a significant Maxwellian compo- 
nent di?/di?R decreases roughly exponentially with in- 
creasing Sr, with the decrease being steepest for small 
m x . In this case, for both the HM and IGEX data, for 
small m x the lowest energy bins, which have the high- 
est numbers of observed events, are most constraining, 
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whilst for larger m x the higher energy bins are the most 
constraining. Care therefore needs to be taken to pro- 
duce overall 90% confidence limits. Using upper confi- 
dence limits on the number of events per bin to produce 
exclusion limits for experiments without background sub- 
traction appears to be more involved than previously re- 
alised, and more sophisticated analysis techniques need 
to be developed. Since the Milky Way halo may well 
be poorly approximated by a standard Maxwellian halo, 
the analysis technique employed should not pre-suppose 
an exponentially declining form for dR/dE^. In fact it 
would be useful, if possible, to formulate an analysis tech- 
nique which decouples the WIMP cross-section from the 
assumed WIMP speed distribution. 
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